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Abstract. The XMM-Newton spectra of the isolated neutron star RXJ0720.4— 3125 obtained over 4.5 years can be described 
by sinusoidal variations in the inferred blackbody temperature, the size of the emitting area and the depth of the absorption 
line with a period of 7.1 ±0.5 years, which we suggest to be the precession period of the neutron star. Precession of a neutron 
star with two hot spots of different temperature and size, probably not located exactly in antipodal positions, may account for 
the variations in the X-ray spectra, changes in the pulsed fraction, shape of the light curve and the phase-lag between soft and 
hard energy bands observed from RX J0720.4— 3125. An independent sinusoidal fit to published and new pulse timing residuals 
from a coherent analysis covering ~12 years yields a consistent period of 7.7±0.6 years supporting the precession model. 
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1. Introduction 

The 8.39 s X-ray pulsar RXJ0720.4 -3 125 was d i scover ed in 
the ROSAT all-sky survey data by iHaberletalJ ill 9971) and 
identified with a faint blue optical star ("iVIotch & Haberfl998^, 
iKulkarni & van Kerkwiik 1998; Kaplan etal. 2 003) which 
show s a proper motion of about 100 mas/yr jMotch et aP 
l2003 >). The bright, soft X-ray source belongs to a small group 
of nearby radio-quiet isolated neutron stars wit h blackbody- 
like thermal X-ray spectra (for recent r eviews see lTreves et alJ 
l2000>. 'Motchl 200 1 1 lHaherll2004 1200.51) . 

Broad absorption lines in the X-ray spectra of most of 
these objects were reported and are usually interpreted as due 
to resonant absorption at the proton cyclotron energy and/or 
bound-bound transitions in H or H-like He ( Haberl et al. 200i 
Ivan Kerkwiik etalJl2004 IZane etalj Eooi iHaberl et al]r2004l 
hereafter H04). This is suggestive of neutron star magnetic 
field strengths of « 10^^-10^'* G, values which are consistent 
with those derived from pul se timing for the period deriva- 
tives of RXJ0720.4-3125 ('Zane et al.' Eool Icropper et alJ 
L2004; Kapl an & van Kerkwiik 200 5a, hereafter K05) and 
RBS1223'( lKapl an & van KerkwiikI l2005bl) if one assumes 
magnetic dipole braking. Moreover, RX J0720.4— 3125 shows 
variations of the absorption line depth with pulse phase (H04) 
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with the line being weakest near intensity maximum and deep- 
est at the declining part of the pulse. 

Among the seven thermal isolated neutron stars discovered 
with ROSAT, RX J0720.4-3125 is unique by showing a grad- 
ual change of the X-ray spectrum on a time scale of years ac- 
companied by an energy-dependent change in the pulse pro- 
file (|de Vries et al. 2004, hereafter D04). In this letter we re- 
port on further XMM-Newton observations of this enigmatic 
source which show that the long-term trends have reversed. We 
present pulse-phase averaged and pulse-phase resolved spec- 
tra from the EPIC-pn instrument. The results of our spectral 
and temporal analysis show strong evidence for a cyclic vari- 
ation of the X-ray properties which strengthens the case of a 
precessing neutron star 

2. XMM-Newton observations 

The soft X-ray source RX J0720.4— 3125 was observed with 
XMM-Newton (Jansenetal. 2001) as calibration target until 
2004 when it was recognised to exhibit a variable spectrum. 
Since then we monitored the spectral evolution of the source. 
Here we utilise the data collected with the EPIC-pn camera 
dStruder et al]l200ll) . The details of the XMM-Newton obser- 
vations with the instrumental setup used for EPIC-pn are sum- 
marised in Table[T] RX J0720.4-3125 is slighdy too bright for 
the full-frame (FF) read-out mode causing photon pile-up ef- 
fects which systematically harden the spectrum (for details of 
this effect see H04). For a relative comparison of the spectra 
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Table 1. XMM-Newton 
RXJ0720.4-3125. 



EPIC-pn observations of 



Orbit 


Observation 


Date 


Setup^ 


Exp. [s] 






0078 


0124100101 


2000 May 13 


FF thin 


44275 




o 


0175 


0132520301 


2000 Nov. 21-22 


FF medium 


22952 


o 




0533 


0156960201 


2002 Nov. 6-7 


FF thin 


25696 


U 




0534 


0156960401 


2002 Nov. 8-9 


FF thin 


27273 






0622 


0158360201 


2003 May 2-3 


SW thick 


51022 




o 
in 


0711 


0161960201 


2003 Oct. 27 


SW thin 


12711 






0711 


0161960201 


2003 Oct. 27-28 


SW medium 


17429 




o 


0815 


0164560501 


2004 May 22-23 


FF thin 


26642 






0986 


0300520201 


2005 April 28 


FF thin 


37101 




1 


1060 


0300520301 


2005 Sep. 23 


FF thin 


34989 






1086 


0311590101 


2005 Nov. 12-13 


FF thin 


33777 







13-05-2000 (rev 0078) 
06-11-2002 (rev 0533) 
08-11-2002 (rev 0534) 

22- 05-2004 (rev 0815) 
28-04-2005 (rev 0986) 

23- 09-2005 (rev 1060) 
12-11-2005 (rev 1086) 




Read-out mode and filter; FF: Full Frame; SW: Small Window. 



we ignore this effect but emphasise that small systematic dif- 
ferences are present between the FF mode spectra and those 
obtained in the faster small window (SW) mode which are free 
of pile-up effects. We do not use the observation with thick fil- 
ter due to the much lower efficiency and different detector re- 
sponse. We screened out strong background flares and list the 
resulting net exposure times (including dead time which is im- 
portant for SW mode) in Table [2 We processed the data using 
the XMM-Newton analysis system SAS6.5 and extracted spec- 
tra and light curves from circular regions with radius of 30". 

2. 1. Pulse-phase averaged spectra 

We extracted and analyzed pulse-phase averaged spectra for the 
10 observations (at different epochs and different instrumen- 
tal setup) following H04 who found that an absorbed black- 
body model with a broad Gaussian absorption line represents 
the spectra of RX J0720.4-3125. We adopted this model and 
fitted the spectra simultaneously with the absorption column 
density (Nh), the line energy and width common to all spec- 
tra. We allowed temperature and the normalizations of black- 
body and line to vary freely between the spectra (except when 
observations were performed shortly after each other as those 
from satellite revolutions 533/534 and the two parts of 711). 
We present the EPIC-pn spectra obtained in FF mode with thin 
filter together with their best-fit model in Fig. [l] which shows 
that the long-term trend of spectral hardening has reversed. We 
report the fit results in Table |2] The common parameters ob- 
tained from the best fit are Nh = 1.01 ±0.03 xlO^o cm"^, a 
line energy of 280±6 eV with a width of cr = 90±5 eV. 

The best fit values for the blackbody temperature (kT), the 
line equivalent width (EW) and the inferred emitting area (ra- 
dius for a circular region with an assumed source distance of 
300 pc) as function of time are shown in Fig. |2l All parame- 
ters can be described with a sinusoidal variation and we derive 
a period of 2580±180 days for the variation in temperature. In 
the fit we do not include the values obtained from the S W mode 
spectra because of the possible systematic difference (pile-up, 
see above). Among the three parameters the temperature can be 
probably best constrained from the X-ray spectra. It is also not 
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Fig. 1. Pulse-phase averaged EPIC-pn spectra of 
RX J0720.4— 3125 from the seven observations with the 
same instrumental setup (FF read-out mode, thin filter). The 
same detector efficiency for all the spectra allows a direct 
comparison and demonstrates the long-term spectral changes. 
The softest spectrum (the uppermost at low energies) was 
obtained in May 2000, while the hardest (the lowest at low 
energies) is that from May 2004. After May 2004 the spectra 
increased monotonically at low energies. 

Table 2. Spectral analysis of pulse-phase averaged X-ray spec- 
tra. 



Orbit 


kT [eV] 


EW [eV] 


Flux'^' [erg 


cm~^] 


0078 


86.6±0.4 


-5.02±4.5 


0.98x10 


-11 


0175 


86.5±0.5 


+8.68±7.7 


1.01x10 


-11 


0533/534 


88.3±0.3 


-21.5±2.6 


0.98x10 


-11 


0711/711 


91.3±0.6 


-73.7±4.9 


1.19x10 


-11 


0815 


93.8±0.4 


-72.4±4.7 


1.02x10 


-11 


0986 


93.5±0.4 


-68.3±5.2 


1.02x10 


-11 


1060 


93.2±0.4 


-67.4±4.3 


1.02x10 


-11 


1086 


92.6±0.4 


-67.5±3.5 


1.04x10 


-11 



Throughout the paper errors are given for a 90% confidence level, 
'^'observed flux is given in the energy band 0.1—2.4 keV. The sys- 
tematically higher flux from satellite revolution 711 is caused by the 
systematic difference of FF and SW mode (FF pile-up losses). The 
statistical uncertainties on the flux are less then 2x10^^* erg cm~'^ 
s^^, which are negligible compared to the stability of the instrument 
of about 1% between observations. 



clear if there are additional variations in the absorption line en- 
ergy and width which we assume to be constant in our current 
analysis. Therefore, we also fit a sine to the variations in EW 
and emission radius, but fix the period at the value found from 
the temperature variation. 

2.2. Pulse-phase variations 

We assigned pulse phases to all detected events using the X-ray 
timing ephemerides inferred by K05 (the "All Data" solution) 
and folded the light curves in different energy bands to produce 
pulse profiles. Examples from two XMM-Newton observations 
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Fig. 2. Variation of temperature, line equivalent width and ra- 
dius of the emitting region of RX J0720.4— 3125 derived from 
a simultaneous fit to the EPIC-pn spectra using an absorbed 
blackbody model with a broad absorption line at 280 eV. The 
sine wave with a period of 7.1 years indicates the best fit to 
the data derived from the FF mode (marked with filled circle) 
observations. For the sine fit of EW and radius we fixed the 
period at the value derived from kT. SW mode data from Oct. 
2003 (open circle) are not used in the fits. 



are presented in Fig. |3] As was found by D04 the pulse pro- 
file became deeper with time. In 2004 the shape of the profile 
changed in particular in the hard band and the hardness ratio 
was higher on average and showed more modulation as com- 
pared to the observation in May 2000. To investigate the spec- 
tral evolution as function of pulse phase we divided the pulse 
into five phase intervals of equal length with the first interval 
from phase 0.0 to 0.2. 

To avoid any systematic shifts due to different instrumental 
setups, we use only the seven observations in FF mode with 
thin filter Similarly to the analysis of phase-averaged spectra, 
we performed a simultaneous fit with the same model to the 
5x7 spectra. Nh, the line energy and width were treated com- 
mon to all spectra and were found to be consistent within the 
errors with the values derived from the phase-averaged spec- 
tra. In Fig.l^the derived line EW is plotted versus temperature 
kT. For each observation the evolution of the two parameters 
during the X-ray pulse follows an ellipse-like track (sampled 
by five points from our finite number of phase intervals) in the 
kT-EW plane. The evolution proceeds counter-clockwise; the 
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Fig. 3. Pulse profiles of RX J0720.4-3 125 in two different en- 
ergy bands (soft S: 0.1 2-0.4 keV; hard H: 0.4-1.0 keV) together 
with the hardness ratio HR=H/S. The top panels are obtained 
from the observation in May 2000 and the bottom panels from 
that of May 2004. 

point marked with a circle indicates the phase interval 0.0-0.2. 
Several remarkable features are seen in Fig.|3 1) The variation 
in kT was smaller during the first observations, consistent with 
the ^2.5 eV value reported by H04. During the later obser- 
vations the amplitude in the kT variation increased to ~6 eV, 
almost as large as the long-term change of ~8 eV seen in the 
phase-averaged spectra (Fig. |3. 2) The amplitude in the line 
EW variation is ^-^40 eV and did not change significantly be- 
tween the observations. 3) The long-term trend reversal of the 
evolution seen from the phase-averaged spectra (Fig. |2j after 
the May 2004 observation is similarly seen at all pulse phases. 
We note that during the pulse variation of the first observation 
the line is formally detected as emission line in phase interval 
0.8-1.0, i.e. before reaching the intensity maximum. Also dur- 
ing the second observation there is an indication that the line 
in the phase-averaged spectrum is in emission rather than ab- 
sorption. However, we can not exclude that this is caused by 
uncertainties in the calibration which may result in a system- 
atic shift of the EW by ~10 eV (for all spectra). 

3. Discussion 

The medium resolution X-ray spectra obtained by the EPIC-pn 
instrument show that the long-term spectral evolution discov- 
ered by D04 has reversed. The spectra became softer following 
the observation of May 2004 when the phase-averaged spec- 
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Fig. 4. Equivalent width of the absorption line vs. temperature 
kT derived from the FF observations with thin filter. Lozenges 
denote the values derived from the phase-averaged spectra. 
During the pulse the parameters evolve counter-clockwise, the 
circle marks phase 0.0-0.2. 



trum was hardest. All spectra can be well modeled by an ab- 
sorbed blackbody continuum with a broad absorption line of 
Gaussian shape. The spectral variations both on time scales of 
years and during the 8.39 s pulse can be explained by vari- 
ations in the inferred blackbody temperature, the depth of the 
absorption line and the size of the emitting area. In the simplest 
approach we did not allow changes in the energy and width of 
the line. If the line is caused by cyclotron absorption both pa- 
rameters may be expected to vary. This requires a more detailed 
investigation of the spectra also involving the higher resolution 
RGS data which is work in progress. 

Blackbody temperature, size of emitting area and absorp- 
tion line depth (equivalent width) derived from the phase- 
averaged spectra follow a sinusoidal variation with a period 
of 7.1 years. Although the XMM-Newton observations do not 
yet cover a full period, it should be noted that the tempera- 
ture of 79±4 eV derived from the ROSAT spectrum (H04), 
obtained in Sept. 1993, ~6.6 years before the first XMM- 
Newton observation is also consistent with a periodic behav- 
ior of RX J0720.4— 3125. The involved time scale is strongly 
suggestive of a scenario with a freely precessing neutron star, 
as first suggested by D04. The observed period of ~7.1 yr is 
then naturally interpreted as the precession period of the neu- 
tron star 

In that case the precession period should also be visible in 
the residual phase shifts seen in the coherent timing analysis 
of the spin evolution of RX J0720.4-3125 (K05). Therefore, 
we fitted the phase residuals as inferred by K05 with a sinu- 
soidal function, including the three EPIC-pn data points for 
the new XMM-Newton observations of 2005. The additional 
data points formally exclude the cubic model used by K05 and 
favour the sinusoidal model. We obtain a preliminary period of 
2830±220 days (7.7±0.6 years) from the phase residual anal- 
ysis. This period is somewhat longer, but consistent with that 
inferred from the spectral analysis. The three ROSAT obser- 



vations extend the time coverage to more than 12 years, i.e. 
more than 1.5 precession cycles. The results of a detailed tim- 
ing analysis which directly includes precession in the spin fre- 
quency model will be presented in a forthcoming paper 

A picture based on thermal emission from two hot spots on 
the surface of a freely precessing neutron stars appears indeed 
promising in explaining many of the peculiar characteristics 
of RX J0720.4— 3125 which have been a challenge so far. In 
this scenario the long-term change in temperature is produced 
by the different (phase-averaged) fractions of the two spots 
which enter into view as the star precesses. In order for such 
a model to work, the two emitting regions need to have dif- 
ferent temperatures and sizes, as it has been recently proposed 
in the case of anothe r member of the same class (RBS1223, 
ISchwope et all 12005'). The nearly, but not exactly sinusoidal 
pulse profile (Fig. |3] D04) and the - to first order - success- 
ful fit of a phase-connected timing solution indicates that the 
poles can not be exactly antipodal, but must be predominantly 
in the same (East- West) hemisphere. 

This is again similar to what was found for RES 1223. 
During the first XMM-Newton observation (May 2000), kT ~ 
85 eV and the cooler, larger spot was predominantly visible, 
while four years later (May 2004), precession brought into 
view mostly the second, hotter and smaller spot, increasing the 
temperature to ~95 eV. This may explain the observed tem- 
perature/emitting area variations and their anti-correlation (see 
Fig. 13 . It should also be noted here that the total flux we see 
from the source stays quite constant, i.e. the two poles con- 
tribute about the same amount. The presence of two different, 
non-antipodal spots may also account for the observed changes 
in the pulsed fraction, light curve shape and for the apparent 
(spin) phase-lag seen between the soft and hard band. 

In Fig. |5] we show the results obtained for such a 
model with a nu merical code adapted from that discussed by 
IZane & TuroUal ll2006ll . First the phase-resolved spectrum is 
computed, having chosen values of the temperature and size 
of the two spots (Ti 2, ^1,2), and of their relative angular 
displacement d^, to match the main observational features of 
RX J0720.4— 3 125. We then compute the phase-averaged spec- 
trum and repeat the calculation at different precession phases 
for a chosen value of the precession angle a. The two angles 
X and ^ denote the inclination of the line-of-sight with respect 
to the precession axis and that between the centre of one spot 
and the rotation axis, respectively. We assume isotropic black- 
body emission from the caps and no attempt has been made 
to model the absorption feature. Relativistic ray-bending is in- 
cluded. Since the present calculation has mainly illustrative 
purposes, we do not account for the detector response and in- 
terstellar absorption nor do we look for parameter fine-tuning. 

As can be seen from Fig. |5] the model is capable of re- 
producing the observed variations of the blackbody tempera- 
ture and emitting area, and also their phase anti-correlation. 
The anti-correlation of the hardness ratio with the total inten- 
sity is also recovered, and the maximum of the hardness ra- 
tio shifts ahead in phase. The light curve modulation evolves 
and the shape changes, becoming more sharply peaked at cer- 
tain epochs. Although no proper fitting of the observed light 
curves and spectra has been attempted at this stage (work in 
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Fig. 5. Left panel: the variation of the blackbody temperature and radius (both normalised) as a function of the precession phase. 
Right panel: the evolution of the light curve (0. 12— 1 keV band) and hardness ratio with precession phase; full, dashed and dash- 
dotted lines refer to phases 0.05, 0.4 and 0.75, respectively. Here a = 20°, Ti = 80 eV, T2 = 100 eV, sin 61 = 0.8, sin 62 = 0.6, 
6*0 = 160°, X = 75° and^ = 5°. 



progress), our preliminary analysis strongly supports the idea 
that RX J0720.4— 3125 is a precessing, nearly aligned rotator, 
seen almost equator-on. First indications seem to favour preces- 
sion angles ^ 10°, larger than those found for radio-pulsars. 
We point out that the modulation is sensitive to the value of 
M/R (2GM/Rc2 = 0.42 in our example of Fig. EJl. Although 
present uncertainties prevent us from reaching any firm con- 
clusion, future improvements in the modelling and new obser- 
vations may allow us to derive constraints on the equation of 
state. 

The variations we see in the absorption line could also be 
related to the different properties of the two emitting spots. Two 
misaligned caps are suggestive of a magnetic field structure 
more complex than a core-centred dipole. A different field con- 
figuration could be then the origin of the different temperature 
and size of the spots. In both RBS1223 and RX J0720.4-3125 
the inferred size of the hotter spot is smaller than that of the 
cooler spot, probably because the former is more confined due 
to a steeper surface gradient of the magnetic field. A further 
possibility is that the line changes (both with spin and pre- 
cession phase) arise because of geometrical effects. The rays 
which reach the observer are at different angles with the mag- 
netic field at different phases. This means that the contribution 
of ordinary and extraordinary photons to the total spectrum is 
phase-dependent. Since extraordinary photons contribute most 
to resonant absorption at the proton cyclotron energy, one ex- 
pects the line properties to change with phase. 

A detailed analysis of the implications of the precession of 
RX J0720.4— 3125 for the physics of neutron star interiors is 
outside the scope of this letter We note, however, that the star 
distortion e = (I3 — Ii)/Ii = Pspin/Pproc ~ 4 x 10~^ is larger 
than that reported for radio-pulsars (e.g. Jon es & Andersso*!] 
I2OOI ■ Akgun et al. 2006.) but smaller than that for HerX-1 
jKetsaris et alJ200oi) . 
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